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Abstract 

We show that minuscule entropic forces, on the order of 100 fN, can prevent the formation 
of DNA loops-a ubiquitous means of regulating the expression of genes. We observe a tenfold 
decrease in the rate of Lacl-mediated DNA loop formation when a tension of 200 fN is applied 
to the substrate DNA, biasing the thermal fluctuations that drive loop formation and breakdown 
events. Conversely, once looped, the DNA-protein complex is insensitive to applied force. Our 
measurements are in excellent agreement with a simple polymer model of loop formation in DNA, 
and show that an anti-parallel topology is the preferred LacI-DNA loop conformation for a generic 
loop- forming construct. 
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Since Jacob and Monod's groundbreaking work on gene regulation, the lac operon has 
become a canonical example of how prokaryotic cells regulate the expression of genes in 
response to changes in environmental conditions [1]. The lac operon is responsible for the 
efficient metabolism of lactose in Escherichia coli bacteria ensuring that enzymes capable 
of digesting lactose are produced only when needed. The lac repressor-mediated DNA loop, 
which is formed when tetrameric lac repressor protein binds to two lac operator sites si- 
multaneously, is an important part of this gene regulatory network and is crucial for the 
repression of lac genes [2j. Long range genetic regulation by DNA looping, however, is not 
unique to the lac operon, but appears in a variety of contexts within prokaryotes, such as 
the ara or gal operons and is ubiquitous within eukaryotes 3] . While the biochemistry of 
these processes is generally well understood, the mechanics of the assemblvand breakdown 
of protein- mediated DNA loops has only recently garnered much attention 4| . In this paper, 
we investigate the role that tension in the substrate DNA plays in the formation and break- 
down of protein-mediated DNA loops, and conclude that loop formation is acutely sensitive 
to entropic forces on the hundred-femtonewton scale. 

Protein-mediated DNA loop formation is driven by thermal fluctuations in the DNA 
which bring distant operators close enough for loop closure by a protein. However, it is quite 
surprising that the magnitude of these fluctuations, which one can estimate as ksT/l p pa 80 
fN where l p = 50 nm is the persistence length, is much smaller than the typical piconewton 
forces that arise in the intracellular environment, from molecular motors or DNA-cytoskeletal 
attachments for example. This observation has led to predictions that forces as small as a 
few hundred femtonewtons are sufficient to reduce the loop formation rate by more than 
two orders of magnitude [5-7|. Given that the cellular environment is thought to regularly 
subject DNA-protein complexes to large static or fluctuating forces, the cell must either use 
mechanical pathways to regulate genetic function, or compensate for the effects of tension 
to ensure the stable control of gene expression. 

To experimentally study the effects of tension on the kinetics of DNA looping, we used 
optical tweezers in conjunction with tethered-particle motion (TPM) measurements to inves- 
tigate the formation and breakdown of Lacl-mediated DNA loops under a constant stretching 
force. We report three main results: First, the rate of loop formation is extremely sensitive 
to applied tension resulting in a tenfold decrease in loop formation when increasing the 
tension from 60 to 183 fN. Second, the lifetime of the looped state appears to be completely 
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FIG. 1: (a) Tethered DNA is trapped in the linear region of an optical potential (dashed line 
indicates laser focus). The figure also represents a kinetic model of DNA looping, (b) Raw experi- 
mental recording of LacI- mediated DNA looping. The looped/unlooped threshold is chosen at the 
minimum between the two state distributions displayed in a histogram of the binned data. 



unaffected by forces as large as 183 fN. Third, our measurements strongly suggest that 
the anti-parallel conformation is the dominant topology of a generic Lacl-mediated DNA 
loop [8] . We mechanically attenuate the thermal fluctuations that drive loop formation and 
breakdown, and measure the associated changes to the looping and unlooping rates, by em- 
ploying a variety of optical tweezers that differs from the more conventional tweezers setup. 
Constant-force axial optical tweezers stretch the molecule away from the surface and trap 
the attached microsphere slightly below the laser focus in an approximately linear region 
of the optical potential. This provides effectively a constant force in the axial direction 
that does not change when the protein binds to or dissociates from the DNA (see Fig. la). 
Details of this set-up are described in Ref.js]. Because of this novel optical tweezers setup, 
we have been able to measure the formation and breakdown rates of Lacl-mediated loops as 
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a function of applied tensile force in the femtonewton range. 

The DNA samples used in this study were prepared in a similar way to that of other 
TPM experiments 9_|. We surface-tethered a 1316-bp ds-DNA molecule with two symmetric 
lac operators spaced 305 bp apart and then attached an 800 nm polystyrene microsphere 
to the other end, which was trapped within the linear regime of the optical potential. The 
total tension in the DNA was carefully calibrated to account for the applied optical force ^ 
and volume exclusion effects arising from entropic interactions between the microsphere and 



the coverslip [lOj. The looping and unlooping lifetimes were measured under four different 
forces: 60 (±5), 78 (±6), 121 (±9), and 183 (±15) fN in the presence of 100 pM of Lad 
protein. In each measurement, the surface-tethered ds-DNA molecule was stretched by a 
contant force while the CCD camera captured defocused images of the tethered microsphere 
at a frame rate of 100 fps. The looped and unlooped states of the DNA molecule, which 
correspond to different axial positions of the microsphere, can be measured by analyzing the 
resulting images, as shown in Fig. lb. By directly observing changes in the axial position 
of the microsphere, the temporal resolution for detecting loop formation and breakdown 
events in our experiment is as low as 300 ms, an order of magnitude better than conven- 
tional TPM. However, as we decrease the applied tension, it becomes increasingly difficult 
to resolve changes in the size of the microsphere, and at zero optical force, we must resort 
to conventional TPM. Moreover, even in the absence of an optical force, a residual entropic 
force from excluded volume effects remains. For this reason, we were not able to obtain a 
direct measure of the force free loop formation and breakdown rates. 

The data was analyzed by first extracting the elapsed times between loop formation and 
breakdown events from time traces like the one in Fig. lb. Then, for each force condition, the 
lifetime of each state was determined from a fit to the cumulative probability distribution, 
as shown in Fig. [21 The resulting distribution displayed by the looped state is well fit by a 
single exponential function 

P{t ] T) = l-e- t '\ (1) 

with time constant r. However, the data of the unlooped state is poorly fit to a single 
exponential function, but is well fit to a biexponential function 

P2(t; n,r 2 ) = cP(t; n) + (1 - c)P(t; r 2 ), (2) 

with time constants T\ and r 2 , and dimensionless fitting parameter c. Results of the fits are 
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FIG. 2: Cumulative probability distributions of the observed durations that the DNA molecule 
remains in the looped and unlooped state under increasing force conditions (left to right). The 
data of (a) the looped state and (b) the unlooped state are fit to the single exponential function 
and a biexponential function respectively. 

shown in Table 1. 

One of the most striking features of the data in Fig. [2] is that the dissociation time 
constant of DNA loops is unaffected by increasing the force from 60 to 183 fN. This result 
is in contrast to the force dependence of the time necessary to form a loop, which increases 
significantly with only a modest increase in applied tension. To interpret these observations 
quantitatively, we begin by applying a kinetic model for the underlying processes of protein 
binding, unbinding, loop formation, and breakdown, as illustrated in Fig. la. This is the 
simplest model of the kinetics that is both consistent with our data and what is currently 
known about LacI mediated looping. 

If we collect all time intervals that start at a loop formation event (L) and end at a loop 
breakdown event (Si), then, within this ensemble, simple first-order kinetics are given by 
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the process 

L ^ S u (3) 

where L is the looped state, Si is the state of the DNA with only one operator bound 
to a protein, and ku is the unlooping rate. Therefore, the time dependent probability of 
unlooping is 

S 1 (t) = l-e- ku \ (4) 

which corresponds to the fit function (Eq.[TJ for the lifetimes of the looped state. The kinetics 
of loop formation, however, are more complicated because there are different unlooped sub- 
states that cannot be distinguished within our experiment. We start by collecting all time 
intervals that begin at an unlooping event and end with the formation of a loop. The kinetics 
may be represented as 

S 2 ^S 1 ^ L, (5) 

fc_ 

where S\ is the state of one vacant and one occupied operator and may directly convert to 
the looped state L at a rate k L , or remain unlooped and convert to state S 2 at a rate 
State S*2, however, is an alternate configuration with both or neither operator occupied by 
a protein, which is not able to directly form a loop, but may convert to state Si at a rate 
k + . With the initial condition Si(0) = 1, the first-order kinetics above may be solved for 
the time- dependent probability of forming a loop 

L{t) = 1 — — — \{k — ku + a)e~ t/Tl - (k - k L - a)e~ t/T2 ] , (6) 
za 

where k = k + + a = [(k + ku) 2 — 4/c+A^] 1 / 2 , and the time constants are defined as 
Ti = 2/ (k + ku — a) and r 2 = 2/(/c + &l + «)• Equation |6] is again a biexponential distribution 
and corresponds to the fit function of Eq. [21 Therefore, we can unambiguously extract the 
four rate constants in our kinetic model. The results are shown in Table 2 and plotted in 
Fig. |3j Our main observation is that, within the uncertainties of our measurements, k+, 
k-, and ku are independent of force whereas ku is acutely force-sensitive on the hundred- 
femtonewton scale. This is consistent with the conventional expectation that the rate of 
conversion between the unlooped states Si and 5*2 does not vary significantly as a function 
of the stretching force on this scale. The insensitivity to applied force of the unlooping rate 
ku can be explained by considering the binding energy of the LacI protein to the DNA, whose 
disassociation from the lac binding site is necessary to break a loop. With a binding energy 
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TABLE I: Fits to the cumulative probability distributions. 



Force (fN) 


60 ± 5 


78 ±6 


121 ±9 


183 ±15 


Looped 


7-00 


20.8 ±0.3 


22.5 ±0.2 


24.1 ±0.3 


24.8 ±0.7 


Unlooped 


ri(s) 


3.0 ±0.2 


4.0 ±0.42 


9.7 ±0.4 


14.5 ±1 


T 2 (s) 


31 ±8 


54 ± 10 


91 ± 10 


101 ±6 


C 


0.77 ±0.03 


0.77 ±0.02 


0.74 ±0.02 


0.38 ±0.03 



of 10 



-19 



11] and an operator region that spans ~ 20 bp, the minimum force needed to 
remove the protein from the operator is ~ 10 pN, which is several orders of magnitude greater 
than the tension we applied. It is clear then why the looped state is relatively insensitive 
to mechanical tension. On the other hand, the sensitivity of the looping rate to such small 
forces is quite striking and potentially rich in implications. Since the characteristic force 
that results from thermal fluctuations of ds-DNA is approximately 80 fN, and since DNA 
looping is a result of thermal fluctuations, femtonewton forces can clearly impact the loop 
formation process. 

Quantitatively useful models of loop formation must explicitly consider the orientation of 
the operators along the DNA in the looped state, as the exact geometry of the loop matters 
significantly. Such theories were developed by Blumberg et al. |6| and, independently, by Yan 
et al. [7J. In this paper, we use the model developed by Blumberg et al. so begin by finding 
the difference in the force dependent contributions to the free energy between a looped and 
a stretched length of DNA: AF = Fi(f, 6) — Fs(f). The kink angle 9 is defined as the angle 
between the tangent vectors of the DNA at the operator sites of the protein-DNA complex. 
A relation for the excess contribution to the free energy as a function of kink angle, imposed 
on the DNA by the loop, is given by: 

F = 4/ 1 / 2 (l-cos(fl/4)) 

L 1 + 12/- 3 /2(i _ cos(6'/4))/(l + cos(tt - 6»))' { ) 

where the free energy is in units of ksT and the force / is in units of the characteristic force 
for thermal fluctuations, f c = k^T jl, p « 80 fN. 

An analytic relation for the free energy of a stretched segment of DNA is given by the 
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difference between the potential energy of a worm-like chain (WLC) and the work done by 
tension: 

hx 2 ( 1 \ 

' / ^ (8) 



4 \{\-xf 

where li is the loop length of the DNA and x is the relative extension of the DNA in units of 
l p . We can now calculate the characteristic time necessary to form a loop under an applied 
force using the principle of detailed balance 

^ = e~ AF , (9) 

where tq is the characteristic time at zero force. 

X-ray studies of the co-crystals of LacI protein bound to short operator fragments have 



revealed the structure of the DNA protein complex 12]. These results impose constraints 
upon, but do not fully determine the topology of the DNA loop. The preferred direction in 
which the DNA enters and leaves the looped complex remains unsettled, but the correspond- 
ing topologies are either anti-parallel conformations, with a kink angle of approximately 150°, 
or parallel conformations, with a kink angle of 30°. 

To fit the data, we calculate Tf = from Eq. [9] as a function of force using the WLC 
model to provide the relative extension x in Eq. [HJ Since we cannot directly measure the 
force free lifetime To, we use this as a single adjustable parameter to fit the curves. The 
value for To is given by a least squares fit to the data. We then generate a curve for both the 
anti-parallel and parallel conformations and see, from Fig. El that the anti-parallel topology 
is more force-sensitive than its parallel counterpart. Our data suggest that the anti-parallel 
conformation is the dominant topology of a generic Lacl-mediated DNA loop. In conclusion, 
our results establish that very small forces, on the order of a hundred femtonewtons, can 
control the assembly of the regulatory protein-DNA complex necessary for expression of 
the lac gene. On the other hand, once formed, the looped complexes are quite stable 
and cannot easily be disrupted by tension in the substrate DNA, giving the system much- 
needed robustness. Thus, it appears more than likely that mechanical pathways can control 
transcription through the application of tiny forces that are generated by other intracellular 
processes. We hope that the development of force measurement techniques inside living cells 
will lead to the identification of such pathways. We also conclude that such regulatory forces 
would likely act on the assembly process of these complexes, but not their breakdown. 
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FIG. 3: Measured values for /c_(A), fc+(o), ki(M), and ku{o) (see Table 2). The looping rate 
ki is fit by the theoretical predictions for the anti-parallel (solid line) and parallel (dashed line) 
topologies illustrated in the insert. 



TABLE II: Rates extracted from kinetic rate equations. 



Force (fN) 


60 ±5 


78 ±6 


121 ±9 


183 ± 15 


Kinetic 
Rates 
(l(r 3 /s) 


ku 


48.1 ±0.6 


44.5 ± 0.4 


41.5 ±0.6 


40 ±1 


k L 


262 ±7 


196 ±6 


79 ±1 


32 ±1 


k_ 


61 ±9 


48 ± 5 


21 ±2 


25 ±4 


k + 


40 ± 10 


24 ± 5 


14 ±2 


21 ±3 
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